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Abstract 

A simple 1-dimensional model is constructed for polymer motion. It exhibits 
the cross-over from reptation to Rouse dynamics, through gradually allowing 
hernia creation and annihilation. The model is treated by the density matrix 
technique which permits an accurate finite-size-scaling analysis of the behavior 
of long polymers. 

1 Introduction 

It is known for some time that the universal properties of long poljTiiers can be obtained 
from stochastic lattice models, which in themselves are rather crude representations 
of the intricate polymeric motion. The reason is that long polymers are critical PP 
and in critical systems the universal properties are independent of the microscopic 
details. There are basically two modes of motion for polymers. One is reptation, 
which is the mechanism for polymers dissolved in a gel and to a lesser extent for dense 
polymer melts. Here the polymer is strongly confined and the main degree of freedom 
is motion inside the confining tube. The other mode applies for dilute solutions where 
the polymers can also move freely sideways. This is usually called Rouse dynamics. 
One can easily envision situations where a mix of these two mechanisms is present. 

The interesting aspect is that the two modes have different and unusual dynamical 
exponents. In polymer motion the dynamic exponent is related to the renewal time 
r. It increases as r ~ A^^, where z is the dynamic exponent and A a measure for the 
length of the chain. Whereas many models have the dynamic exponent z = 1, showing 
isotropy between time and space, reptation has an exponent z = 3. This value has 
been a bit controversial, since viscosity measurements point at z = 3.4, whereas the 
theories agree on z = 3. The discrepancy has recently been removed by an accurate 
finite-size-scaling analysis using the density-matrix technique (DMRG) introduced by 
White pij. Rouse dynamics on the other hand is found to have the exponent z = 2. 
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It is therefore interesting to study the cross-over between these two mechanisms and 
to find out how a mixing-in of Rouse dynamics changes the dynamical exponent from 
z = 3 to z = 2. 

The most convenient model for reptation is the Rubinstein-Duke (RD) model de- 
signed by Rubinstein |lj and extended by Duke [3] , by introducing a driving field. The 
mobile units, the reptons, only move along the tube that the chain traces out in the 
lattice. The advantage of the model is that the dimension d of the lattice, in which 
the polymer chain is embedded, becomes a parameter, which infiuences the behavior 
of the ends of the chain but not of the bulk. As this parameter d is one of the details, 
having no infiuence on the universal properties, one often studies [0] the case d = 1, 
although the RD model becomes somewhat artificial in a 1-dimensional embedding. 

There are two forms of sideways motion of the reptons. When a cell is occupied by 
three reptons the middle one can enter a neighboring cell without crossing a barrier. 
This is called hernia creation and the opposite process is hernia annihilation. The 
other forms of sideway motion imply that the chain crosses a barrier. These are the 
typical motions allowed in Rouse dynamics. Within the spirit of the physics of the RD 
model, hernia creation and annihilation should be allowed, but that makes the model 
essentially more difficult. For instance, the role of the embedding dimension cannot be 
simply reduced to a parameter d, infiuencing only the ends of the chain. 

The usual argument to omit the hernia creation and annihilation, is that these 
processes do not alter the universal properties. This is likely to be true in larger d, 
where hernias become a fraction of the possibilities for the chain, but in ci = 1 they 
are of major importance as we will show in this paper. In fact the hernia creation 
and annihilation mimic the role of Rouse dynamics in a 1-dimensional embedding and 
therefore it is a convenient mechanism to study the cross-over from reptation to Rouse 
dynamics. 

The dynamic exponent z is obtained from the gap in the spectrum of the Master 
Equation. Apart from this gap another interesting quantity is the diffusion coefficient. 
We obtain this from the model by studying the drift velocity in the limit of a weak 
driving field. For reptation the diffusion coefficient decays as N~'^ for chains of length 
A^, while for Rouse dynamics the diffusion is speeded up to A^^^. Next to the cross- 
over of the dynamic exponent, we study in this paper the cross-over of the diffusion 
exponent. 

2 The model 

The model is a 1-dimensional chain of A^+1 reptons, connected by A^ links, {yi, ■ ■ ■ , y^). 
The links are either in the forward direction, = 1, or in the backward direction 
yi = —1, or have the value y^ = 0. The cases yi = ±1 are considered as taut links, 
while ?/j = is a slack link or an element of stored length. The basic motion rule is the 
hopping of this stored length unit along the chain, by interchanging with taut links. 
If it moves in the forward direction, its transition rate is biased by a factor B > 1, 
while the hopping rate in the backward direction is decreased by the factor B^^ < 1. 
The biases represent an external field driving the reptons of the chain. At the end of 
the chain the links may change from slack to taut and vice versa, thereby adding or 
subtracting an element of stored length, again with a bias depending on the direction 
of the transport of length. These motion rules form the much studied Rubinstein-Duke 
(RD) model. Our new element is that we allow a neighboring pair of opposite taut links 
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to change into a pair of slack links and vice versa. We describe this as the annihilation 
viz. creation of a hernia. The transition rate for hernia creation/annihilation is h, 
multiplied with a bias based on the sign of motion of the middle repton of the hernia. 

Without hernia motion the RD model is a typical model for reptation. The tube, 
which is the sequence of taut links, can only be changed from the ends. This is a slow 
process, since the taut links in the bulk have to wait till they happen to drift to one of 
the ends, before they can change their value. Simple counting tells that the inner taut 
links need at least N'^ repton moves, if they could renew themselves in a systematic 
way. The change of configuration is however a diffusive process in configuration space 
and therefore the average renewal time is measured in single repton moves, or A^^ 
in chain updates. So the reptation renewal time r ~ A^^. Obviously hernia creation 
and annihilation speed up the renewal of the chain and the point of this note is to see 
how they can overtake the reptation mechanism. 

A similar global argument ^ yields that the pure RD model (without hernia cre- 
ation/annihilation) has a drift velocity decaying as A^~^, leading to an asymptotic A^~^ 
behavior for the diffusion coefficient. 



3 The Master Equation 

Our model is, as all the hopping models, governed by the Master Equation for the prob- 
ability distribution P(Y) where Y stands for the complete configuration ■ ■ ■ ,?/Ar). 
It has the form 

^^^^ = Y: [W^(Y| Y')P(Y', t) - W{Y'\Y)P(Y, t)] ^ M(Y, Y')P(Y', t). (1) 

The Ws are the transitions rates and the matrix M contains the gain terms (in the off- 
diagonal elements) and the loss terms (on the diagonal). Conservation of probability 
implies that the sum over the columns of the matrix vanishes. So the matrix has a zero 
eigenvalue and the eigenfunction corresponding to this eigenvalue is the stationary state 
of the system, to which every other initial state ultimately decays. The matrix is non- 
symmetric, due to the bias, which gives different rates to a process and its inverse. Thus 
one has to distinguish between left and right eigenf unctions. The left eigenfunction 
belonging to the zero eigenvalue is trivial (all components equal); the problem is to 
find the right eigenfunction as the stationary state probability distribution. 

The renewal time is given by the slowest decaying eigenstate. Thus the gap in the 
spectrum near is the inverse renewal time. All eigenvalues must have of course a neg- 
ative real part, otherwise probability would grow unlimited. The form stresses the 
similarity to quantum mechanical problems. Indeed the linear structure of the poly- 
mer chain makes it a 1-dimensional quantum problem, however, with a non-hermitian 
hamiltonian. Our approach to the solution exploits this analogy by applying the so- 
called DMRG method for quantum problems to find the properties of the transition 
matrix M. In previous publications Tj the application of this method to polymer 
motion has been described in detail. Here we present only the results. 

We confine ourselves to the renewal time and the diffusion coefficient. The renewal 
time is usually defined at zero driving field. Also the standard diffusion coefficient 
refers to zero driving field . However, to determine the diffusion coefficient we must 
turn on an infinitesimal field and compute the drift velocity. This can be done by 
expanding the Master Equation with respect to the field. The field enters in the bias 
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B which we represent as 

S = exp(e/2), (2) 

where e is a dimcnsionless parameter measuring the field strength. Then we expand 
the Master Equation in powers of e 

M = Mo + eMi + ---, P(Y) = Po(Y) + ePi(Y) + --- (3) 

and obtain the equations 

MoPo = 0, MoPi = -MiPo- (4) 

The first equation is trivially fulfilled by a constant Pq{Y), since the matrix Mq is 
symmetric and the right eigenvector becomes equal to the trivial left eigenvector. The 
second equation is a set of homogeneous linear equations for the components of Pi(Y). 
It is soluble, since the right hand side of the equation is perpendicular to the left 
eigenvalue (which remains true to all orders is e) . So we can make the solution definite 
by requiring that it is also orthogonal to the trivial left eigenvector. Pi(Y) yields the 
lowest order drift velocity and the diffusion constant follows by the Einstein relation 
as 

D--(^] . (5) 




Figure 1: The renewal time as function of the length of the chain for various values of 
the hernia creation/annihilation rate h. 
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Figure 2: The diffusion exponent a; at as function of tfie lengtfi of tfie cliain for various 
values of the hernia creation/annihilation rate h. 

4 Scaling exponents 

One of the advantages of the DMRG method is that it calculates the properties, e.g. 
the gap, for a growing length of the chain. In principle the method allows to go to 
any length, but it is in practice limited by instabilities and computational time. We 
have speeded up the process by using the field inversion symmetry in linear order in e, 
both for the gap and the diffusion coefficient. This make the results very well suited 
for a finite-size-scaling analysis. We convert the gap as function of to a renewal time 
t{N). In Fig. [T]we present the local exponent zn, defined as 



The DMRG method gives values, accurate enough, such that the small differences 
in (jni) do not spoil the accuracy. The various curves correspond to different values of 
h. Previously we found that it is most suggestive to plot r as function of N~^^'^ . This 
applies indeed for the case = 0, but for non-zero values of a plot agains A^~^ gives 
more straight curves. Some features are noteworthy: 

• Chains of the order of ~ 100 are not yet in the asymptotic regime [7j. So 
there are large corrections to scaling. This is the origin of the earlier mentioned 
controversy between theory and experiment. In particular the plateau in the 
/i = curve (the pure reptation case) may easily lead to the conclusion that the 
exponent has settled on the (too large) value. 



Zn — 



lnr(iV + 1) - lnr(A^ " 1) _ d\nT 
ln(A^+ 1) -ln(iV- 1) ~ d\nN 



(6) 
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• The influence of small values of h is quite strong for long chains in particular for 
small values of h. We come back on this point when we discuss the cross-over 
behavior. 

• The asymptotic behavior of the exponent (for — > oo) differs for h = from 
all the other curves. While the theoretical value Zqo = 3 for reptation, is quite 
compatible with the data, it is definitely excluded for the curves h ^ 0. They 
clearly point to the common value Zoo = 2, which is characteristic for Rouse 
dynamics. 

In Fig. 121 we plot in the same way the local exponent xn for the diffusion coefficient, 
defined as 

~ ln(A^ + 1) - ln(A^ - 1) ~ ~ d\nN' ^ ' 

The picture has a similar message as the previous one. It is clear that, without hernia 
motion {h = 0), the exponent evolves towards the value 2, while for any non-zero 
value of h, it aims at the value 1. Again one has large corrections to scaling. These 
corrections make it impossible to determine the exponent from In-ln plots. Only due 
to the high accuracy of the DMRG method one can derive exponents from formulae 
like ® or (0). 

In Fig. 121 we have made a plot of ln(r/A^^) and — ln{DN). In both cases the 
asymptotic values — > oo are plotted as function of In h. As one sees the curves are 
fairly straight, with a slope -0.55, in the domain where the data are most accurate. For 
very small values of h we see in the renewal data a somewhat smaller slope, a trend 
which is also detectable in the diffusion data on closer inspection. 
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Figure 3: In- In plots of the renewal time and the diffusion coefficient as function of h 



6 



5 Crossover scaling 



The point of crossover scaling is to represent the data for various values of h in one 
single curve. Anticipating the asymptotic values of the two regimes: h —>■ and a fixed 
h ^ 0, the following representation is adequate for the renewal time. 

T{N,h) = N^g{h^N). (8) 

The connection with the previous representation runs via the relation 

rflnr _ d\ng{h^N) 
dlnN ~ ^ ^ d\n{hm) ■ 

We expect the cross-over function g{x) to be expandable for small arguments as 

9ix)=go + giX^ (10) 

and for large arguments as 




Inserting the asymptotic behavior (fTT|) into ^ we obtain 

\n{T/N^) = \ng^i-e\nh + ---, (12) 

where the dots refer to corrections of order 1/A^. So the slope in Fig. El gives the value 
of 6. The value g{0) can be derived from a plot of tN~'^ versus N~^. We find the value 
g[0) ~ 0.2. In Fig. lUwe have plotted the scaling function as a function of h^'^^N. 
The observed data collapse is the proof for cross-over scaling. The deviation for small 
argument in Fig. 0] are due to short chains. The scaling curve should aim, for small 
arguments, at the value g{0) ~ 0.2, which we deduced from the h = curve. 
In Fig. El we plot similarly the diffusion coefficient in the form 

D{N, h) = N-^f{h^N) (13) 

with the same value 6 = 0.55. As one sees the collapse is excellent. The cross-over 
scaling function / approaches again a finite value at x = 0. In view of the data 
for the diffusion coefficient at /i = we have /(O) = 0.4, which is quite consistent 
with the behavior of the h > curves. For large arguments, f{x) should behave as 
/(x oo) ~ X. 

We hesitate to claim that the cross-over scaling exponent differs from the value 
6 = 1/2, which certainly gives a less perfect data collapse. An argument in favor 
of = 1/2 is based on the simple estimate of the times to remove a hernia for the 
two mechanisms. As we mentioned, pure reptation requires N"^ single repton moves 
to refresh the chain as a whole. On that timescale the hernias (in total of order N) 
in the chain are annihilated and replaced by others. So it takes repton moves 
to forget a hernia by reptation. On the other hand, direct change of a hernia by 
creation or annihilation goes with a rate h/N. The fastest process dominates and the 
competion is controlled by the ratio of the rates {h/N) / (l/N^) = hN"^. So the crossover 
scaling function should be a function of the ratio hN^. It might well be that the real 
asymptotic value for ^ = 1/2 and that we see in the window, where we have data, an 
effective exponent. This looks similar to the story of the renewal exponent itself, which 
also was estimated as 3.4, while the true theoretical value is 3. As we mentioned we 
see in Fig. Ola tendency to a smaller slope for the very small h, which supports this 
possibility. 
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Figure 4: The cross-over function g as defined in (jH)) as function of the argument h^-^^N . 



6 Discussion 

We have presented a simple model which demonstrates the cross-over from reptation to 
Rouse dynamics. In the Rubinstein-Duke model the links in the direction of the field 
and those against the field cannot interchange and this makes reptation a slow process. 
In our 1-dimensional model, hernia annihilation and creation, allow the two types of 
links to interchange and therefore these obstacles can be overcome. In that sense they 
play the same role as the tube changes which are typical for Rouse dynamics. 

In a paper by Sartoni and Van Leeuwen, [Hj the 1-dimensional reptation with hernia 
creation and annihilation, has been connected to a simpler model of two types of 
particles, which move independently of each other along the chain. They also conclude 
that the diffusion coefficient decays as iV~^, but they have to stick to a hernia creation 
and annihilation rate equal to the hopping rate of the reptons. Here we could vary 
this rate at will and therefore study the cross-over behavior. In a forthcoming paper 
we have related their findings to the recently introduced necklace model CD] . 
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